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Abstract 

The elastic and large-angle quasi-elastic scattering reactions were studied with the same nucleus- 
nucleus potential proposed for describing fusion reactions. The elastic scattering angle distributions 
of some reactions are reasonably well reproduced by the proposed Woods-Saxon potential with 
fixed parameters at energies much higher than the Coulomb barrier. With an empirical barrier 
distribution based on the modified Woods-Saxon potential and taking into account the influence of 
nucleons transfer, the calculated quasi-elastic scattering cross sections of a series of reactions are 
in good agreement with the experimental data. 
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I. INTRODUCTION 



Heavy-ion quasi-elastic scattering and fusion reactions at energies around the Coulomb 
barrier have been extensively studied in recent decades, since they provide an ideal oppor- 
tunity to obtain the information of nuclear structure and nucleus-nucleus interaction and 
to explore the mechanism of heavy-ion reactions at near barrier energies which is of great 
importance for the synthesis of super-heavy nuclei . Based on the 

quantum tunneling concept, it is thought that the quasi-elastic scattering (a sum of elastic, 
inelastic scattering, and transfer channels) is a good counterpart of the fusion reaction in 
the sense that the former is related to the reflectionprobability of a potential barrier while 
the latter is related to the penetration probability [2J. In addition, it has been shown that 
the fusion barrier distribution generated by the coupling of the relative motion of the nuclei 
to internal degrees of freedom can be extracted from precisely measured fusion excitation 
functions jj, 4]. The similarity of the barrier distribution can be extracted from large-angle 
quasi-elastic scattering excitation functions [5| which can be more easily measured than the 
fusion excitation functions |lo|. Therefore, it is expected that both the fusion and quasi- 
elastic scattering cross sections of a heavy-ion reaction at energies around the Coulomb 
barrier can be unifiedly described by the same nucleus-nucleus potential. However, in re- 



cently published papers [ll|, [12[, Mukherjee et al. found that the Woods-Saxon nuclear 
potential can not simultaneously reproduce precise fusion and elastic scattering measure- 
ments of 12 C+ 208 Pb Hj, and Zamrun and Hagino found that the depth parameter of the 
Woods-Saxon potential for describing the fusion cross sections of 16 0+ 144 Sm have to be read- 
justed to reproduce the experimental quasi-elastic scattering cross sections of this reaction 
with the same coupled- channels framework [3] . To solve this discrepancy, it is necessary to 
find a nucleus-nucleus potential for a unified description of the scattering and fusion data 
in heavy-ion reactions. In addition, for giving satisfying predictions of quasi-elastic scatter- 
ing cross sections for unmeasured reaction systems, it is required to find a nucleus-nucleus 
potential to describe quasi-elastic scattering reactions systematically. 

Studies of quasi-elastic scattering reactions and transfer processes, especially of the be- 
havior of the transfer probabilities as functions of the distance of closest approach or the 
incident energies have attracted a lot of attention. Some investigations show that the 
semiclassical method is suitable for describing the heavy-ion scattering at large reaction 
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distance 14, Q, Q. The transfer probability is expressed as an exponential function 
of the distance between the reaction partners P tT oc exp(— 2aR c ) [l3| in the semiclassical 
approximation, where a is the transfer form factor and R c is the distance of closest ap- 
proach between two nuclei. The exponential dependence on R c is a characteristic property 
of tunneling 15) . At energies below the barrier the experimental slopes are generally in 
good agreement with the predictions of the model for one-nucleon transfer. At higher en- 
ergies, the measured slopes deviate from the calculated values, which is often referred to 
as "slope anomaly" in addition to other types of anomalies which are in connection with 
slopes obtained in two particle transfer reactions [nj]. Some experiments show an energy 
dependence of the slopes, and a clear trend of a decrease of slope paramet function of 

nn 

increasing energy was found in |16l. 1171] . The transfer probabilities at below barrier energies 
have been extensively studied while a theoretical model for describing the slope parameters 
and the transfer probabilities at energies near and above the Coulomb barrier has not been 
well established yet. The study of the transfer probability in the latter energy region is still 
required. In addition, it is interesting to explore the relation between the transfer probabili- 
ties and the barrier distribution since the transfer probabilities generally peak in the vicinity 
of the barrier energies [lo| . 

In [3| we proposed a modified Woods-Saxon potential model based on the Skyrme energy- 
density functional together with the extended Thomas-Fermi approach. This model was first 
proposed in [9] and a large number of fusion reactions have been described satisfactorily well 
with an empirical barrier distribution which is based on the calculated entrance channel 
potential. In this work, we try to describe the heavy-ion elastic and quasi-elastic scattering 
with the same potential for describing the fusion reactions. The paper is organized as follows: 
In Sec. II, the theoretical model for the description of the elastic and quasi-elastic scattering 
is introduced. In addition, some calculated results are compared with experimental data. 
The summary and discussion are given in Sec. III. 
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II. THEORETICAL MODEL FOR ELASTIC, QUASI-ELASTIC 
SCATTERING AND FUSION 



In this section, we first briefly introduce the modified Woods-Saxon potential and the 
elastic scattering is studied with the potential. In addition, the empirical barrier distribu- 
tion is briefly introduced for describing fusion reactions. Then, the quasi-elastic scattering 
and the transfer probabilities are described with the empirical barrier distribution. Some 
calculated results are also presented in this section. 

A. Modified Woods-Saxon Potential and Elastic Scattering at above Barrier 

Energies 



In [18| we proposed a Woods-Saxon potential model based on the Skyrme energy-density 
functional together with the extended Thomas-Fermi approach. The nucleus-nucleus inter- 
action potential reads as: 

V(R) = V N (R) + V C (R). (1) 

Here, Vn and Vq are the nuclear and Coulomb interactions, respectively. We take Vc(R) = 
e 2 ZiZ 2 /R, and the nuclear interaction Vn- 

Vn{R) = l + eM(R-Ro)/aY (2) 



with 11 



A 1/3 A 1/3 

V = u [1 + k(I 1 + I 2 )} / 3 2 (3) 

A 1 + A 2 



and 

R = roiA 1 / 3 + A] /3 ) + c. (4) 

Ii = (Ni — Z\)jA\ and l 2 = (N% — Z 2 )/v4 2 in Eq.(3) are the isospin asymmetries of the 
projectile and target nuclei, respectively. In this potential, the depth of the potential V 
depends on the reaction system and the isospin asymmetries. To distinguish it from the tra- 
ditional Woods-Saxon potential (with three parameters) in which the depth of the potential 
is independent of the reaction system and the isospin asymmetries, we call the proposed 
potential (with five parameters) "modified" Woods-Saxon potential. The parameters of the 



TABLE I: Parameters of the modified Woods-Saxon potential, 
ro (fm) c (fm) uo (MeV) k a (fm) 
1.27 -1.37 -44.16 -0.40 0.75 



modified Woods-Saxon (MWS) potential 18] are determined by the entrance channel po- 
tentials of 66996 reactions obtained with the Skyrme energy-density approach and are listed 
in Table 1. 

The proposed nucleus-nucleus potential is based on the frozen density approximation. The 
time-dependent Hartree-Forck (TDHF) calculations show that the nucleus-nucleus potential 
depends on the incident energy at energies close to the Coulomb barrier and when the 
center-of-mass energy is much higher than the Coulomb barrier energy, potentials deduced 
with the microscopic theory identify with the frozen density approximation 20j. We test 
the modified Woods-Saxon potential for the description of heavy-ion elastic scattering at 
energies much higher than the Coulomb barrier, since the reaction time is relatively short 
and the frozen density approximation seems to be applicable at these energies. Based on 
the optical model, we solve the Schroedinger equation for a given nucleus-nucleus potential 
using the traditional Numerov method to obtain the partial-wave scattering matrix that is 
used to describe the elastic scattering data [2l| . The real and imaginary parts of the optical 
potential adopted in the calculations are described by the modified Woods-Saxon potential. 

We have calculated the elastic scattering angular distributions for the reactions 
i2 C+2 08 Pb) i6 O+ 208 Pb) i2 C+ 90 Zr and i6 0+ 63 Cu at different laboratory energies. The cal- 



culated results (so 
(squares) 2^, 25, 



id curves) are shown in Fig. 1, and the corresponding experimental data 



26 



27 



28j are also presented for comparison. The experimental data 



of the four reactions at different energies are reasonably well reproduced by the modified 
Woods-Saxon potential in which the potential parameters are fixed. 

We further test the MWS potential for the description of heavy-ion fusion at above barrier 
energies. At these energies, the fusion cross section is usually described by the classical 
formula 



(?fus(E c . m .) = tiR) (1 - B/E, 



(5) 



with the fusion radius Rf and the height of the fusion barrier B. Fig. 2 shows the fusion 
excitation function of 16 O+ 208 Pb. Taking B to be the barrier height B (78.72 MeV) of the 
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FIG. 1: (Color online) Elastic scattering angular distributions for the reactions 12 C+ 208 Pb, 
16 O+ 208 Pb, 12 C+ 90 Zr and 16 0+ 63 Cu at different laboratory energies. The solid curves and the 



squares denote the calculated results with the modified Woods-Saxon potentia 



tal data, respectively. The experimental data are taken from Refs. 



24, 
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26, 



and t 



27 



le experimen- 



28]. 



modified Woods-Saxon potential, the fusion cross sections at above barrier energies can not 
be reproduced by the Eq.(5) (see the dash-dotted curve in Fig. 2). In order to describe the 
fusion cross sections satisfactorily, we introduced an empirical barrier distribution to take 
into account the mult i- dimensional character of a realistic barrier due to the coupling to 
internal degrees of freedom of the binary system in our previous paper {9]. We proposed an 
effective weight function for describing the barrier distribution, 
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FIG. 2: (Color online) Fusion excitation function of the reaction 16 O+ 208 Pb. The dash-dotted 
and solid curves denote the calculated results with Eg .(5 ) by taking B = Bq and B = -B m . p ., 



29(] . The sub- figure shows the effective 



respectively. The solid circles denote the experimental data 
weight function of the reaction. 

where D avT (B) = (Di(B) + D 2 (B))/2 and B x is the left cross point of D X {B) and D 2 (B). 
Di(B) and D 2 (B) are two Gaussian functions jl], [lsj which depend on the barrier height B 
of the modified Woods-Saxon potential. The effective weight function D e ^ of the reaction 
16 O+ 208 Pb is shown in the sub-figure of Fig.2. Taking B to be the most probable barrier 
height -B m .p. (74.43 MeV) according to the -D e fr, the fusion cross sections at above barrier 
energies are reproduced reasonably well (see the solid curve in Fig.2). With the empirical 
barrier distribution, the fusion cross sections and the mean barrier heights of a large number 



of reactions can be reproduced well 



0,Q, 



22|. 



From the above discussion, one finds that for the heavy-ion elastic scattering at above 
barrier energies the modified Woods-Saxon potential which is based on the frozen density 
approximation gives nice results. But the fusion cross section of the same reaction system can 
not be described well with the potential and the barrier distribution needs to be introduced 
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to reproduce the fusion data. In recently published paper [23|, the authors proposed two 
optical potentials for describing the reactions 12 C+ 208 Pb and 16 O+ 208 Pb, respectively. Both 
the elastic scattering and fusion data can be satisfactorily described with the potentials at 
energies around the Coulomb barrier. At energies much higher than the Coulomb barrier (for 
example, E^ = 192 MeV for 16 O+ 208 Pb), the elastic scattering data can not be reproduced 
well by the potential. In this work, we aim to find a nucleus-nucleus potential for describing 
the reactions systematically. 

B. Description of Large-angle Quasi-elastic Scattering 

As a good counterpart of the fusion reaction, the large-angle quasi-elastic scattering is 
studied to explore the nucleus-nucleus potential. In this work, we explore the influence 
of the empirical barrier distribution proposed for the fusion reactions on the large-angle 
quasi-elastic scattering. 

It is thought that the quasi-elastic differential cross section can be expressed as a weighted 
sum of the eigenchannel elastic differential cross sections under the adiabatic and iso- 



centrifugal approximation 30|, l3jj . Similar to the description of fusion with the empirical 
barrier distribution, we describe the large-angle quasi-elastic scattering cross section with 
the effective weight function D c r(B) at energies around the Coulomb barrier, 



(E c .m.) — PeS + -Pcorr, (7) 



da R 
with 

1 poo J 

Pe S = Tr D cS (B)^(E c . m .,B)dB, (8) 
F Jo da R 

and P C orr is a small correction term. is the ratio of the elastic cross section a e \ to the 
Rutherford cross section ctr. F is a normalization constant F = J D c ^(B)dB. Within 
the semi-classicalperturbation theory, a semi-classical formula for the backward scattering 
(9 = tt) is given |d, lad], 



Where the nuclear potential V^(R C ) is evaluated at the Coulomb turning point, 

v N(Rc) M-^f\ or^-^j ) , do) 

V Rf J \l + exp[(R c - R )/a\J 



with the distance of the closest approach between two nuclei R c = Z\Z 2 e 2 / E c , m .. a is the 
diffuseness parameter of the nuclear potential. Z\, Z 2 , E c _ m , denote the charge numbers of 
the projectile and target nuclei and the center-of-mass energy, respectively. Rf and hw are 
the barrier position and curvature of the modified Woods-Saxon potential, respectively. 

The correction term P corr in Eq.(7) takes into account some effects in the quasi-elastic 
scattering that are not involved in the empirical barrier distribution (which was proposed for 
describing fusion reactions). In this work, we assume that the correction term mainly comes 
from nucleons transfer. In principle, the transfer process also affects the fusion process and 
the effect of nucleons transfer may have been implicitly taken into account in the empirical 
barrier distribution. However the influence of nucleons transfer on the quasi-elastic scattering 
may differ from the influence on the fusion process, and a small correction term seems to be 
required. 

For the quasi-elastic scattering, we first investigate the dependence of the transfer 
probabilities on the incident energies. The transfer probabilities Ptr can be written as 
P tr = (da tT /dQ)/(da^/dQ), where da tT /dfl is the transfer cross section [141 ] . At energies 
below the barrier the transfer probability increases with increasing incident energies be- 
cause the distance of closest approach becomes small, leading to an increase in the nuclear 
overlap. At above barrier energies, on the contrary, the transfer probability decreases with 
increasing energies since the increased overlap results in more dissipative collisions which 
finally results in fusion [16]. For the transfer at energies below the barrier, we describe the 
transfer probability using the traditional semi-classical method which is mentioned previ- 
ously in the introduction section. In this work, only the one-neutron transfer channels are 
taken into account for simplicity in the calculation of the transfer probability at sub-barrier 
energies. For the transfer at energies above the barrier, we find the transfer probabilities 
are close to the Gaussian function D 2 of the empirical barrier distribution. For example, in 
Fig. 3 the measured transfer probabilities for the reaction 16 0+ 232 Th are compared with 
the corresponding Gaussian function D 2 of this reaction at energies above the barrier. The 
experimental data are in good agreement with D 2 at above barrier energies. 
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FIG. 3: (Color online) Transfer probability of 16 0+ 232 Th. The squares denote the measured 
transfer probabilities [16] including the channels of lp, lpln, la, 2p and 2pln transfers to the target 
nuclei. The crossed curve denotes the Gaussian function D 2 in the empirical barrier distribution. 

Based on the above discussion, we write the transfer probability: 



Ptr{E c .m.) — / 



Pq exp(— 2aR c ) 

D2(E c .m.) 



E c , m , — B n 



> B h 



(11) 



B m < E r rn < Bh 



with the strength factor f = 1 MeV. Pq is a normalization constant given by Pq 



D 2 (Bq)/ exp(— 2aZiZ 2 e 2 / B ) with the transfer form factor a = ^j2[iEb/h 2 . /x is the re- 
duced mass of the transferred nucleons, Eb is the effective binding energy of the transferred 



nucleons 



331 ] . Bq is the barrier height of the modified Woods-Saxon potential, B m is the 



(12) 



mean barrier height of the barrier distribution, 

/ B D e ff(B) dB 
* m jD cS (B)dB ■ 

We take B^ = 2Bq — B m in this work. In order to have a smooth function for the transfer 
probability P tr from sub-barrier energies to above barrier energies, we introduce a function for 
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P t r in the energy region B m < E c _ m _ < with a factor F = D ^ — y exp(— 2a ZiZ 2 e 2 / B m ) 
to link the two functions P exp(— 2aR c ) and D 2 describing the P tI at sub-barrier and at 
above barrier energies, respectively. 

We assume P corr ~ P tr . Both the fusion and quasi-elastic scattering cross sections of 
a series of reactions have been studied with the proposed approach in this work. Fig. 
4 to Fig. 6 show the calculated quasi-elastic scattering and fusion cross sections for the 

Sm, 16 0+ 92 Zr, 16 0+ 186 W, 32 S+ 208 Pb and 16 0+ 116 Sn. The 



reactions lb O+ 144 Sm, lb O 



experimental data 



34 



154 



38 



39 



40] are also presented for comparison. The solid circles 



and squares denote the measured fusion cross sections <jf us and large-angle quasi-elastic 
scattering cross sections, respectively. The solid curves in (a) and (c) of Fig. 4 to Fig. 6 
denote the calculated results for <jf us with the proposed empirical barrier distribution (see 
details in Refs.jsl, 1^]). The crossed curves in (b) and (d) denote the calculated quasi- 
elastic scattering cross sections with Eq.(7). The dashed curves denote the results for P e g, 
i.e. the contribution of the empirical barrier distribution to the quasi-elastic scattering. 
We find that both the fusion excitation functions and the quasi-elastic scattering excitation 
functions of the six reactions can be satisfactorily well reproduced. In Fig. 7 we compare the 



measured quasi-e 
35^ i6 0+ 232 Th 



astic scattering excitation functions (squares) of the reactions 12 C+ 142 Nd 



16], 16 0+ 64 Zn |3J| and 32 S+ 110 Pd [37] with the calculated results with 



Eq.(7) (crossed curves). The calculated quasi-elastic scattering cross sections of the four 
reactions 12 C+ 142 Nd, 16 0+ 232 Th, 16 0+ 64 Zn and 32 S+ 110 Pd are in good agreement with 
the experimental data. Fig. 4 to Fig. 7 indicate that the modified Woods-Saxon potential 
together with the empirical barrier distribution can simultaneously describe the quasi-elastic 
scattering and fusion of a number of reactions reasonably well. 



III. CONCLUSION AND DISCUSSION 

In this work, we have studied the heavy-ion elastic and large-angle quasi-elastic scatter- 
ing with the same nucleus-nucleus potential proposed for describing fusion reactions. The 
elastic scattering angle distributions of a series of reactions at energies much higher than the 
Coulomb barrier can be reasonably well reproduced by the modified Woods-Saxon potential 
which is based on the frozen density approximation systematically. With the same potential 
the fusion cross sections at above barrier energies can not be reproduced and an empirical 
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FIG. 4: (Color online) Fusion cross sections and quasi-elastic scattering cross sections as a function 
of energy for the reactions 16 0+ 144 Sm and 16 0+ Sm. The solid circles and squares denote the 
measure fusion cross sections <7f us and quasi-elastic scattering cross sections, respectively. The solid 
curves in (a) and (c) denote the calculated results for cjf us . The crossed curves in (b) and (d) denote 
the calculated results with Eq.(7). The dashed curves denote the results for P e R . The experimental 
data of fusion and quasi-elastic scattering are taken from Ref.[J| and [5], respectively. 

barrier distribution (which takes into account the coupling of other degrees of freedom) is 
required to reproduce the fusion data. It seems that the coupling of other degrees of free- 
dom to the relative motion of the nuclei is obvious in heavy-ion fusion processes whereas the 
frozen density approximation is applicable in the elastic scattering process at energies much 
higher than the Coulomb barrier. With the empirical barrier distribution function based on 
the modified Woods-Saxon potential, the fusion cross sections of a series of reactions have 
been well reproduced. Further, with the same empirical barrier distribution and taking into 
account the correction term that mainly comes from the nucleons transfer, the calculated 
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FIG. 5: (Color online) The same as Fi^ 4, but for the reactions 16 0+ 92 Zr and 16 0+ 186 W. The 



experimental data are taken from Refs. 



large-angle quasi-elastic scattering cross sections of these reactions are in good agreement 
with the experimental data. 
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experimental data are taken from Refs. 



. 4, but for the reactions 32 S+ 208 Pb and 16 0+ 116 Sn. The 
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